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^ ■ Abstract 

o" 

^^ i We present an approach for the calculation of spin density distributions for molecules 

2 \ that require very large active spaces for a qualitatively correct description of their elec- 

tronic structure. Our approach is based on the density-matrix renormalization group 
f*-^ . (DMRG) algorithm to calculate the spin density matrix elements as basic quantity for 

CN I the spatially resolved spin density distribution. The spin density matrix elements are 

directly determined from the second-quantized elementary operators optimized by the 
F^ . DMRG algorithm. As an analytic convergence criterion for the spin density distri- 

bution, we employ our recently developed sampling-reconstruction scheme [J. Chem. 
Phys. 2011, 134, 224101] to build an accurate complete-active-space configuration- 
QJ ■ interaction (CASCI) wave function from the optimized matrix product states. The 

spin density matrix elements can then also be determined as an expectation value em- 
^ ■ ploying the reconstructed wave function expansion. Furthermore, the explicit recon- 

O ■ struction of a CASCI-type wave function provides insights into chemically interesting 

^ ! features of the molecule under study such as the distribution of a- and /3-electrons in 

terms of Slater determinants, CI coefficients, and natural orbitals. The methodology 
O i' is applied to an iron nitrosyl complex which we have identified as a challenging system 

for standard approaches [J. Chem. Theory Comput. 2011, 7, 2740]. 
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1 Introduction 

In quantum chemistry, the electronic structure of molecules is described by either ab initio 
wave- function methods or density- functional theory (DFT). For large molecular systems such 
as transition metal complexes, however, wave-function based methods are rarely employed 
due to the corresponding high computational cost (for counterexamples, see Refs.'^'^. Hence, 
the application of DFT became instrumental in theoretical studies of mechanisms in metal- 
mediated catalysis.'SKIll Yet, the treatment of open-shell system^^ISlIISl g^^^^j (near-) degenerate 
states remains a challenge for DFT.^ Failures of approximate exchange-correlation density 
functional in predicting properties of open-shell systems have been traced to the delocal- 
ization error and static correlation error^'^ which are rooted in an inappropriate behavior 
of the energy with respect to fractional charges and fractional spins .'^ In addition to the 
difficult prediction of ground states from states of different spin,'I^'22H30] gp^^ density distri- 
butions considerably depend on the approximate exchange-correlation density functional if 
transition metal complexes containing noninnocent ligands are considered. ^^ Qualitatively 
correct spin density distributions are difficult to obtain within the standard Kohn-Sham 
formalism that has not been formulated to also produce accurate spin densities.'^ 

However, accurate spin densities are desirable for various reasons. (1) In electron paramag- 
netic resonance spectroscopy (EPR),'^the spin density is the central quantity on which EPR 
parameters explicitly depend. ^^ Obviously, reliable spin density distributions are important 
for an accurate calculation of EPR properties, but this remains a difficult task to achieve 
for theoretical chemistry.'^SHSSl (^2) The question which approximate exchange-correlation 
density functional yields sufficiently accurate spin densities remains inconclusive.^^ ^° If ac- 
curate reference spin density distributions were available, a more detailed analysis of the spin 
density distribution in terms of spin density difference plots could be used as a qualitative 
and quantitative benchmark for the validation of approximate exchange-correlation density 
functionals. (3) According to the Hohenberg-Kohn theorem,'^the spin density is not needed 
to calculate the electronic energy or any other expectation value. However, in open-shell sys- 
tems it is often introduced as an additional variable which leads to a spin-DFT formalismpl 
first introduced by von Barth and Hedin.^^ In spin-DFT, the spin density becomes a fun- 
damental quantity and reliable reference spin densities could be used to construct proper 
approximations to the exact exchange-correlation density functional. 

For accurate spin densities in cases for which a DFT description fails, ab initio electron 
correlation methods need to be applied. Pierloot et al presented complete-active-space 
self-consistent-field (CASSCF) studies for large transition metal complexes which provided 
deeper insights into the quality of DFT spin density distributions .'^"Sl xhe large molecular 
size of these systems requires large active orbital spaces, but the standard CASSCF ap- 
proach restricts their dimension which represents the most crucial approximation in such 
calculations. ^° It is therefore important to understand whether the spin density is converged 
with respect to the dimension of the active orbital spaces used so far. This is a task that is 
difficult to study within a standard CAS-type approach. 

In general, up to about 18 electrons correlated in 18 spatial orbitals are computationally 
feasible for standard CASSCF. These limitations may restrict the accurate description of 



electronic structures which could be approved only by enlarging the dimension of the active 
orbital space. Reliable reference spin density distributions for complicated open-shell struc- 
tures as found, for instance, in iron complexes with non-innocent ligands require capabilities 
beyond those of standard correlation methods. 

A different approach for the calculation of correlated ab initio spin densities for large 
molecules was recently presented by Kossmann and Neese^^ who discussed the performance 
of orbital-optimized M0ller-Plesset perturbation theory in calculating hyperfine coupling 
constants for atoms and small molecules. In this approach, isotropic hyperfine constants of 
coupled-cluster singles-doubles quality could be obtained which could be further improved 
by applying spin-component scaling. 

Here, we pursue a diflFerent route for the calculation of ab initio spin densities by applying 
the density-matrix renormalization group (DMRG) algorithm. With the DMRG algorithm, 
introduced by White^l^in 1992, much larger active orbital spaces can be considered beyond 
the limit of, say, 18 electrons correlated in 18 molecular orbitals. It was shown that DMRG 
is capable of providing accurate wave functions and energies even for complicated electronic 
structures (see Refs.^^^^ for reviews). Moreover, we first showed that the DMRG algorithm 
yields reliable relative electronic energies between diflFerent spin states or isomers of transition 
metal complexes and clusters for which DMRG was not meant to work and which are a very 
challenging task for any other multi-reference quantum chemical method^^ (see also Ref.^^ for 
latest results and further references). We shall demonstrate in this work that also accurate 
DMRG spin density distributions can be determined for very large active orbital spaces. 

Recently, we presented a convergence analysis of the spin density distribution for a small 
iron nitrosyl model complex [Fe(NO)]^^ in a field of point charges, which demonstrated that 
medium-sized active orbital spaces are sufficient for quantitatively correct spin densities ."^^ 
However, a quantitative analysis that can explore truly large active spaces is still lacking for 
this complex, which shall therefore be the target system in this work. In such cases, DMRG 
spin densities can be considered as reliable references which can serve as benchmark results 
for approximate exchange-correlation density functionals. 

This work is organized as follows. In section [21 we discuss the spin density matrix and 
its spatially resolved counterpart, the spin density distribution, employing the formalism of 
second quantization. Then, we continue with the introduction of DMRG spin densities. In 
section 12. 2[ we present our approach of approximating the DMRG spin density distribution 
by one from a complete-active-space configuration-interaction(CASCI)-type wave function 
which allows us to compare DMRG spin densities from calculations with different DMRG 
parameter sets. In order to validate our approach, we study the spin densities of a medium- 
sized active orbital space in section [3l This is then extended by considering up to 29 active 
orbitals in section [H Finally, a summary and concluding remarks are given in section [51 



2 Spatially resolved, non-relativistic spin densities 

Since DMRG is based on the second quantized formalism, we briefly discuss how the spin 
density in spatial coordinates can be written in second quantization. In first quantization, 
the operator for the spin density reads 

N 

S^^'''{r) = J2^(^-^^)^^,^^ (1) 

i=l 

where s^^i is the z component of the one-electron spin operator, r^ is the spatial coordinate of 
electron i, and A^ is the total number of electrons in the system. Applying an orbital basis, 
the corresponding operator expression in second quantization is given by 
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= Y,4>;{r)4>,m„ (2) 



p,q 



where p^q run over the total orbital basis {0^} with (j)i{r) representing the spatial part 
of a spin orbital. The operators a-^ and a^^r are the creation and annihilation operators, 
respectively, for an electron of spin a in orbital i. In Eq. ([2]), the spin density operator 
(5^P^^(r) is defined in terms of the spin tensor excitation operators 

Tpq = ^{ala^qa " ^/^^^Z^) (3) 

in the orbital basis (see Ref.'^ for details). The spatially resolved spin density p^P^^(r) is 
calculated as the expectation value of (5^P^^(r), 

p«P-(r) = (^M|<5^P'°(r)|*M) 

= E0;(r)0,(r)(*M|r,,|*M), (4) 

pq 

where |^m) represents some normalized reference state 

|n) = \n1n2 . . .rik) is an occupation number vector with elements n^ G {0, 1}. {n} repre- 
sents the set of aU occupation number vectors constructed from k one-particle states. The 
expectation value on the right hand side of Eq. (|H) is a spin density matrix element Tj 
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T^f^ = {^M\fpg\^M) = ^{"^Mlal^aga - aJ^a^^l^M). (6) 



2.1 DMRG spin densities from second-quantized elementary op- 
erators 

If the reference state |^m) is a DMRG wave function in Eq. (|6j), the corresponding DMRG 
spin density matrix elements Tpq are obtained. The matrix representations of the 

creation and annihilation operators are available in every step of the DMRG algorithm and 
each spin density matrix element can thus be easily determined. 

The operator a^^r^gcr in its matrix representation is calculated as a tensor product for which 
we have to distinguish two different cases. The molecular orbitals p and q are defined either 
(i) on the same or (ii) on different subsystems of the DMRG partitioning of the active orbital 
space into the active (sub)system, its environment (the complementary subsystem), and one 
or two explicitly treated orbitals in between. While the former case is straightforward to 
handle, for an operator expression in the latter case, however, we need to build operators for 
the superblock where all three subsystems, i.e., the active subsystem, the exactly represented 
sites and the environment, are combined as tensor products. 

To illustrate this concept, let us consider two operators ai and a2 defined on three different 
subspaces J^i, ^2, and Ts- Then, the combined subspace T is defined as J^ = J^i (8) ^2 ® ^3, 
where ^ as well as J^i, ^2 and ^s are all subspaces of the A^-particle Fock space ^n. For 
instance, the operator expressions for the combined subspace are given by 

af : ai ® 1^^ ® 1-^^ (7) 

af : A^^ (8)a2 1^3, (8) 

where Aj^. is the anticommutation matrix of the corresponding subspace Ti. For the product 
of two operators we obtain 

af . af = (ai • A^J ® (a2) ® (1^3) , (9) 

where we have used the mixed-product property for the right-hand side of the above equation 
which mixes the ordinary matrix product with the tensor product. All remaining operator 
products can be derived in a similar way. After the spin density matrix is determined, 
the spatially resolved spin density distribution can be calculated from Eq. (|4j). If the wave 
function is real, the spin density matrix is symmetric and the calculation can be speed up 
by calculating the upper triangular part of the spin density matrix only. 



2.2 Spin density from a reduced dimensional CASCI-type wave 
function 

Since CI vectors are in general spars#S^^2^' — if contributions below a predefined threshold are 
neglected — CASCI-type wave functions can be efficiently and accurately projected onto a 
smaller set of Slater determinants, which only represent the most important contributions to 
the wave function expansion. We recently reported the sampling-reconstruction algorithm 
for CASCI-type wave functions defined in a complete active orbital space from a previously 



optimized DMRG wave function (SRCAS algorithm). ^^ An approximate CASCI-type ex- 
pansion I^m) for any wave function |^m) consisting of k one-particle states can thus be 
written as 

I^M) = E^w|fi)' (10) 

{n} 

where the sum runs over all occupation number vectors n living in the sampled subspace 
of the total many-particle Hilbert space. Using Eq. (|6j), we can calculate the spin density 
matrix by substituting the reference state |^m) with the approximate state |^m), 



^^(M[SRCAS]) ^ (^^If^^l^^) = ^{^M\ala,^ - al^a,^\^M) 

= 9 5Z ^w*^{S(fil4-^^--4/3^^/5l^)- (11) 



{n},{m} 

Since the occupation number vectors are orthonormal to each other, the expectation value 
on the right hand side of Eq. (ITT]) can be easily evaluated and we obtain 

7.(M[SRCAS]) _ 1 V^ MM)*MM) f ^ ^ _ ^ X ^ \ f^o\ 

{n},{m} 

where n^ represents the occupation number vector where orbital p lacks one electron with a 
spin. Furthermore, we introduced a phase factor e^r to account for the annihilation operations 
of a^^ acting on the bra-state and a^^r acting on the ket-state. 

Based on this approximate expression for the spin density matrix we can determine spin 
density distributions for subspaces of the many-particle Hilbert space of different dimensions 
and study the sensitivity of the spin density distribution to the number of active-system 
states in DMRG calculations. 



2.3 Measures for spin density comparisons 

For various reasons, we need suitable measures to assess the similarity of different spin 
densities. For instance, such a measure would be required to assess the accuracy of a given 
spin density compared to a reference spin density. 

Monitoring the evolution of the spin density for an increasing number of active-system 
states^^ can illustrate the convergence behavior of the spin density distribution with re- 
spect to the number of active-system states m. Isosurface plots of the difference in spin 
density distributions for calculations with different m- values can only serve as a qualitative 
convergence measure. As quantitative measures, however, we introduce two distances which 
quantify how far two spin densities are apart from each other. Both distance measures are 
defined with the absolute error in the spin density difference distribution. The accumulated 
absolute error Aabs is given by 



^abs 



|pr(r)-pr(r)|dr, (13) 



and the root-square error Ars reads 




spin 



(r) -pr (r)|2dr, (14) 

where p^^^^(r) refers to the spin density distribution corresponding to some calculation i, 
e.g., to a CASSCF or DMRG spin density when different chemical methods are compared, 
or to some parameter sets if different spin densities are determined with the same method. 
If two spin densities p^^^^(r) and p^-^^^(r) are similar, both Aabs and Ars approach zero. For 
accurate ab initio spin densities, we shall require both error measures to be smaller than 
0.005 (Aabs) or 0.001 (Ars), respectively (in view of the results discussed in section IHTTI) . 

A different similarity measure can be applied by employing directly the knowledge of the 
reconstructed CASCI-type wave function expansion. This procedure relies on the closeness 
measure of two quantum states, namely the quantum fidelity.^^^^- The importance and 
potential application of the quantum fidelity within the DMRG framework was first discussed 
by some of u^in the context of quantum error correction and was also utilized in our SRC AS 
approach.l^ Two CASCI-type wave function expansions reconstructed for different numbers 
of DMRG active-system states, rrii and m2, can be explicitly compared by calculating their 
quantum fidelity 

Fm.,m, = \{^^M'^\^^M'^)\' (15) 

as an overlap measure. 



3 A noninnocent model system 

In a previous study, we reported DFT and CASSCF spin density distributions in iron nitrosyl 
complexes as well as for the [Fe(NO)]^+ molecule embedded in a square-planar field of point- 
charges to emulate the one-electron states of the full complexes.— Since DFT spin densities 
of iron nitrosyl complexes remain ambiguous, we choose the small [Fe(NO)]^^ molecule in 
its doublet state for our analysis here. The point charges facilitate a dynamic change of 
the character of the electronic wave function by shortening the distances rfpc of the point 
charges to the metal center. Depending on this distance dpc, both single-reference and multi- 
reference situations can be created for [Fe(NO)]^^. When the four point charges are located 
at a distance of 1.131 A from the iron atom, the electronic structure of the [Fe(NO)]^^ 
molecule represents a single-reference problem, while for d^c = 0.598 A a multi-reference 
case is generated. 

The [Fe(NO)]^+ structure features a Fe-N bond length of 1.707 A and a N-0 bond distance 
of 1.177 A with a Fe-N-0 angle of 146°. The four negative point charges of —0.5 e each 
are located as depicted in Figure [l](e) . Due to the small size of the [Fe(NO)]^+ molecule, 
we can efficiently study the dependence of the spin density distribution on different DMRG 
parameter sets such as the number of DMRG active-system states m. Thereby, we are able 
to define appropriate convergence measures for the spin density in order to reach a predefined 
accuracy. 



3.1 The single-reference case 



As already discussed in great detail in Ref.,HBltlie minimal active orbital space for [Fe(NO)]^+ 
with rfpc = 1.131 A comprises seven electrons correlated in seven orbitals for qualitatively 
reliable spin density distributions. It consists of Fe 3d- {dxy^ dyz^ dxz^ dx2-y2 and ^^2) and 
both NO 7r*-orbitals. As orbital basis in our DMRG calculations, the natural orbitals from a 
CAS(7,7)SCF calculation performed with the Molpro program packagd^ using Dunning's 
cc-pVTZ basis set for all atom^l^ were taken. The one-electron and two-electron integrals 
in the natural orbital basis were also calculated with the Molpro program package. ^^ All 
DMRG calculations reported in this section [3] were carried out with the Zurich DMRG 
program.^ Random noise was added to the density matrix in order to force the mixing of 
configurations that would have not been captured otherwise if the number of active-system 
states m is too small.'^'^ 



We performed DMRG calculations for diflFerent numbers of DMRG active-system states m 
abbreviated by DMRG(x,^)[m] where x corresponds to the number of active electrons and 
y is the number of active orbitals for m renormalized active-system states. Starting with 
m = 16, m is further increased to 32 and 48 until the CAS(7,7)SCF reference energy is 
reproduced for m = 64 active-system states (see Table [I]). Note that the number of active- 
system states needed to reproduce the CASSCF result is very small in this case. This can 
be explained employing concepts of quantum information theory in section [H The DMRG 
calculations reported in this section do not employ these concepts to enforce better conver- 
gence. This decision is deliberately taken in order to produce nonconverged low-m results to 
compare with the m = 64 calculation. We should note that this artifact could be cured by 
the dynamical block state selection (DBSS) procedure,'^'^ while the strong dependence on 
small m- values and the convergence to local minima can be overcome by applying the config- 
uration interaction-based dynamically extended active space (CI-DEAS) procedure.*^ The 

Table I: Ground state energy for [FeNO]^+ surrounded by four point charges at two different 
distance sets rfpc in Hartree atomic units for CAS(7,7)SCF and DMRG(7,7)[m] calculations 
for different numbers of DMRG active-system states m. 



dpc = 1.131 A 


dpc = 0.598 A 


Method 


E/Rartree 


Method 


£;/Hartree 


HF 

CAS(7,7)SCF 


-1392.844 043 
-1392.887 247 


HF 

CAS(7,7)SCF 


-1396.821 220 
-1396.858 313 


DMRG(7,7)[16 
DMRG(7,7)[32 
DMRG(7,7)[48 
DMRG(7,7)[64 


-1392.881 067 
-1392.885 462 
-1392.886 893 
-1392.887 247 


DMRG(7,7)[16 
DMRG(7,7)[32 

DMRG(7,7)[48 
DMRG(7,7)[64 


-1396.762 709 
-1396.818 651 
-1396.840 018 
-1396.858 313 



spin density distributions for our four DMRG calculations (m = 16,32,48,64) are shown in 
Figure [I](a) and were determined as discussed in section 12. 1[ To emphasize the dependence 
on m, the corresponding spin density difference plots with respect to the CAS(7,7)SCF ref- 
erence spin density distribution are displayed. Note that all isosurface plots are shown for 



the same isosurface value of 0.0003. All DMRG calculation yield qualitatively similar spin 
density distribution, only minor quantitative differences can be observed. The CAS(7,7)SCF 
reference spin density can be perfectly well reproduced for m = 64 DMRG active-system 
states and is, hence, not shown in Figure [T](a) . To calculate approximate spin density dis- 
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Figure 1: (a) and (c): Spin density difference plots for DMRG(7,7)[m] spin densities 
calculated for a different number of DMRG active-system states m with respect to the 
CAS(7,7)SCF reference spin densities shown in (b) and (d), respectively, for [Fe(NO)]^^ 
in a quadratic-planar point-charge field. Two different distances are considered, namely 
dpc = 1.131 A, (a) and (b), and d^^ = 0.598 A, (c) and (d). For both distances, the spin 
density could be perfectly reproduced in a DMRG(7,7)[64] calculation and is therefore not 
depicted here. An isosurface value of 0.0003 in (a), (b) and 0.003 in (c), (d), respectively, is 
chosen, (e): Structure of [Fe(NO)]^^ with the four point charges of —0.5 e. 

tributions from reconstructed CASCI-type wave functions, we first have to sample the most 
important configurations of the A^-particle Hilbert space. For this purpose, we applied our 
SRCAS method. ^^ Due to the small size of the active space, the A^-particle Hilbert space 
is spanned by only 1225 Slater determinants and all corresponding CI coefficients can be 
determined directly from the CASSCF reference calculation. In general, similar CI coeffi- 
cients are obtained for all DMRG calculations and the CASSCF reference, i.e., similar wave 
functions are converged resulting in small differences in the spin density distributions. The 
distribution of the CI coefficients is depicted in Figure 1 of the Supporting Information. 



Spin density distributions determined for different sampled subspaces of the A^-particle 



Hilbert are in good agreement with the corresponding DMRG spin density. Note that the 
sampled subspaces are defined by the threshold value of the completeness measure (COM) 
introduced in Ref.^^ with COM = (1 — ^/C|), where / runs over all sampled configura- 
tions with CI coefficients Cj. In general, threshold values of 0.01 to 0.001 turned out to 
be sufficient for obtaining quantitatively reliable spin densities in this single-reference case. 
The corresponding isosurface plots and excitation histograms with respect to the COM are 
summarized in the Supporting Information. 

The spin density difference plots in Figure [l](a) illustrate the convergence of the spin density 
distribution with respect to the number of DMRG active-system states m. The absolute 
error Aabs and the root-square error Ars of the spin density difference distributions provide 
a quantitive measure for the accuracy (see Table [TTI) • The differences in the spin densities 
calculated for m = 48 DMRG active-system states is small compared to the CAS(7,7)SCF 
reference. For 48 active-system states upwards, both Aabs and Ars are below their threshold 
values given in section [Z3l The set of quantum fidelity measures F^-^m^+i for our four DMRG 
calculations with m^ G {16,32,48,64} is {0.980000,0.994395,0.999012}. Increasing m from 
48 to 64 DMRG active-system states corresponds to ^48,64 = 0.999012 which illustrates the 
similarity of both DMRG wave functions and results in reliable spin density distributions 
for m > 48. 

Table II: The absolute error Aabs and the root-square error Ars of the DMRG(7,7)[m] spin 
densities with respect to the CAS(7,7)SCF reference for [FeNO]^+ surrounded by four point 
charges at two different distance sets rfpc employing different numbers of DMRG active- 
system states m. 



Method 


dpc = 


1.131 A 


dpc = 


3.598 A 




Aabs 


Ars 


Aabs 


Ars 


DMRG(7,7)[16] 

DMRG(7,7)[32] 
DMRG(7,7)[48] 
DMRG(7,7)[64] 


0.007678 
0.004392 
0.001397 
1.34-10-^ 


0.002147 
0.001285 
0.000412 
5.40-10-*5 


0.213543 
0.221198 
0.081631 
9.69-10-*5 


0.052168 
0.052144 
0.020418 
3.66-10-'^ 



3.2 The multi-reference case 



A multi-reference character of the [Fe(NO)]^^ molecule can be induced by decreasing the 
distances of the point charges to the iron atom. In the squeezed model complex, the point 
charges are placed at a distance of rfpc = 0.598 A from the iron center in the same configu- 
ration as before. Similar to the single-reference problem, the minimum active orbital space 
considered here comprises seven electrons correlated in seven orbitals. Yet, it consists of four 
Fe dxy^ dyz^ dxz and (i^2 (the d^'^-y^ is excluded due to the compressed point charge environ- 
ment), two NO TT* and one NO a orbital which interacts with the Fe rf^2 orbital. Again, the 
natural orbitals from a CAS(7,7)SCF calculation were taken as orbital basis in our DMRG 
calculations and determined with the MOLPRO program package^ using Dunning's cc-pVTZ 



10 



basis set for all atoms .'^'^ The calculation of the one-electron and two-electron integrals 
in this natural orbital basis was also performed with the MOLPRO program package.'^ All 
DMRG calculations were carried out with the Zurich DMRG program. ^^ As before, we 
performed DMRG calculations for four different numbers of DMRG active-system states m. 
Starting with m = 16, m is further increased to 32 and 48 until the CAS(7,7)SCF reference 
energy is obtained for m = 64 active-system states (see Table [I]). The small-m calculations 
are designed not to reproduce the CAS(7,7)SCF reference for this analysis. Note, however, 
that a small number of active-system states was needed to reproduce the CASSCF result as 
observed in the single-reference problem. 

In Figure [T](d), the CAS(7,7)SCF spin density distribution is shown which is taken as the 
reference distribution, while Figure [T](c) illustrates the spatially resolved differences in the 
DMRG(7,7)[m] and CAS(7,7)SCF spin density distributions. Note that the same isosurface 
value of 0.003 was chosen for all spin densities shown. For small m-values, qualitatively 
different spin density distributions are obtained. The /3-electron density around the nitrosyl 
ligand is underestimated and a dumbbell-shaped /3-electron density is obtained in contrast to 
the cylindric shape of the reference /3-electron density. The a-electron density around the Fe 
atom is underestimated. Increasing m to 48 results in a cylindric /3-electron density around 
the NO-ligand which differs only little from the reference spin density. The spin density can 
be exactly reproduced for m = 64 active-system state for which also the CAS(7,7)SCF refer- 
ence energy is obtained. The convergence properties of the DMRG(7,7)[m] spin density with 
respect to m can be quantified by the Aabs- and Arg-values where significantly large values 
(> 0.005 and > 0.001, respectively) are obtained for spin density distributions determined 
in small-m calculations (see Table [11]) • 

In Figure [21 the distribution of CI coefficients for the DMRG and CASSCF wave functions is 
shown. Since only the position of the point charges has been modified, the A^-particle Hilbert 
space remains spanned by 1225 Slater determinants and all corresponding CI coefficients can 
be determined directly from the CASSCF reference calculation as in the single-reference case. 
Similar CI coefficients are obtained for the DMRG(7,7)[64] calculation and the CAS(7,7)SCF 
reference, i.e., similar wave functions are converged. However, significantly different CI 
coefficients are obtained — as expected — for smaller m-values. In particular, the deviations 
are most significant for configurations corresponding to the largest CI weights. Additional 
information on the distribution of CI coefficients can be found in the Supporting information. 

Similarly to the single-reference problem discussed above, reliable spin densities obtained 
from reduced dimensional CASCI-type wave function expansions can be determined for a 
COM > 0.001 independent of m. A complete collection of spin density distributions for 
different CASCI-type wave function expansions and DMRG parameter sets can be found in 
the Supporting Information. Figure [3] shows the ratio of Slater determinants with respect to 
the complete A^-particle Hilbert space which have been picked up in the sampling procedure 
and sorted by their corresponding CI weights for the DMRG(7,7)[32] and DMRG(7,7)[64] 
calculation. For COM > 0.001, the reconstructed CASCI-type wave function contains the 
major part of the important Slater determinants, while for a further decreased threshold 
value of 10~^ almost all significant Slater determinant have been picked up. Note that the 
sampling procedure was restricted to accept only configurations with CI coefficients larger 
than the threshold value for COM. Although all possible excitations are included in the 
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Figure 2: Distribution of the absolute value of the CI coefficients corresponding to the Slater 
determinants in the DMRG(7,7)[m] calculations with different renormalized active-system 
states m and in the CAS(7,7)SCF reference calculation for [Fe(NO)]^^ surrounded by four 
point charges at a distance of d^c = 0.598 A from the iron atom. All Slater determinants 
are ordered according to the CI weights of the CAS(7,7)SCF calculation. 

CASCI-type wave function in the limit of COM -^ (see also Figure 4 of the Supporting 
Information), the pattern of the CI coefficients of the DMRG(7,7)[64] calculation is different 
from the CI pattern of the CAS(7,7)SCF reference. While large CI coefficients (> 0.0001) 
are reproduced within sufficient accuracy, smaller CI weights are underestimated. The max- 
imum of the curve is shifted towards smaller CI weights < 10~^. Hence, the DMRG algo- 
rithm disregards an exact weighting of unimportant configurations with small CI coefficients 
which is a feature of matrix product and tensor network states where large CI coefficients 
should be reproduced and unimportant configurations are neglect ed."^"^ A complete col- 
lection of excitation histograms for different CASCI-type wave functions can be found in 
the Supporting Information. To quantify the differences in the underlying wave functions 
for our three DMRG calculations employing rrii G {32,48,64} active-system states, we cal- 
culated the quantum fidelity F^-^m^+i which forms in this case a set of overlap measures of 
{0.831887,0.897445,0.955669}. 

We conclude that reliable spin density distributions can be calculated either from converged 
DMRG ground state wave functions or from the reconstructed CASCI-type wave function 
expansions. In particular, a fully converged DMRG wave function is not mandatory to 
obtain qualitatively correct spin density distributions if the CI weights of the most important 
configurations are well reproduced for a given 7n-value. This holds for both the single- 
reference and the multi-reference case. A representative set of Slater determinants, i.e.. 
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Figure 3: CI histogram of the absolute values of the CI coefficients for the Slater deter- 
minants for reconstructed CASCI-type wave function expansions from the DMRG(7,7)[m] 
calculations with different renormalized active-system states m for the [Fe(NO)]^^ molecule 
surrounded by four point charges at a distance of rfpc = 0.598 A from the iron atom. The 
CAS(7,7)SCF reference calculation is also shown for comparison, thr corresponds to the 
threshold value of COM in the sampling-reconstruction procedure and denotes the accuracy 
of the reconstructed CASCI-type wave function. All Slater determinants with CI coefficients 
in an interval as indicated on the abscissa are grouped together. 



the most important ones {\Cj\ 
distribution. 



> 0.001), is sufficient for a qualitatively correct spin density 



4 Spin density distributions for large active spaces 

While we have studied the convergence features of DMRG calculations for small active 
spaces, for which we could obtain an exact CASSCF reference result, we shall now proceed 
to explore territory with DMRG that is not accessible to the CASSCF approach. In our 
recent analysis of CASSCF spin densities for the [Fe(NO)]^+ molecule, ^° the spin density 
distribution was qualitatively converged with respect to the dimension of the active orbital 
space. For quantitatively accurate spin densities we need to increase the dimension of the 
active orbital space so that important iron and ligand orbitals which are missing in the 
standard CASSCF calculations, e.g., the Fe dx2_y2 double-shell orbital, could also be included 
in the active orbital space. Here, we extend the convergence series presented in Ref.^^l by 
considering active orbital spaces containing up to 29 active orbitals. Starting with an active 
orbital space comprising 13 active electrons correlated in 20 active orbitals, the number of 
active orbitals is further increased to 24 and 29, respectively. The two largest active orbital 
spaces do also contain the fifth rf3,2_^2-double-shell orbital which could not be included in all 
CASSCF calculations presented in Ref.^S The [Fe(NO)]^+ molecular structure features the 
same bond distances and angles as presented in Section [3l The four point charges of —0.5 
e are located at a distance of 1.131 A from the metal center in order to properly model the 



13 



square-planar ligand field of the fuU-fiedged complexes in a doublet spin state. 

For all DMRG calculations, the natural orbitals from a CAS(11,14)SCF calculation are 
employed as orbital basis.^^^Hlll Similarly, the CASSCF calculation as well as the calculations 
of the one-electron and two-electron integrals in the natural orbital basis were performed with 
the MOLPRO program package^^ using Dunning's cc-pVTZ basis set for all atoms,-^-^ while 
the DMRG calculations are performed with the Budapest DMRG program.!^ In addition, 
the DMRG orbital orderings were optimized for all three active orbital spaces and the CI- 
DEAS starting guess was performed. Figure [4] displays the corresponding single orbital 
entropies given by 

^(^)^ = - J^cJc.,ilncJc.,i, (16) 

a 

and mutual information determined by 

/,,,• = . (2),,,- -.(1), -.(I),-, (17) 

where i = 1 . . . fc is the orbital index and runs over all k one-particle states, (jUa.i is the 
a eigenvalue of the reduced density matrix of orbital z,^^ while ^(2)^^ is the two-orbital 
entropy between a pair i^j of sites introduced by Rissler et al. to the quantum chemical 
DMRG algorithm."^ Note that the mutual information and single orbital entropies are 
confined to the first ten natural orbitals for all considered dimensions of the active orbital 
space. These natural orbitals are highly entangled and represent the most important orbitals 
comprised in the active orbital space. Therefore, accurate DMRG spin densities can be 
obtained already for a reasonably small number of active-system states. Similar entropy 
profiles can be obtained for smaller dimensions of the active orbital space. 

The number of DMRG active-system states m was set to 128, 256, 512, 1024 and 2048, 
respectively. The ground state energies for all DMRG calculations are summarized in 
Table [TTTl Considering the DMRG(13,20)[m] calculations, an energy convergence of 0.135 
mH (0.4 kJ/mol) is reached with respect to m. For the largest active orbital space, the 
DMRG(13,29)[1024] energy is converged to 1.195 mH (3.1 kJ/mol) when compared to the 
DMRG(13,29)[2048] reference. 



4.1 Convergence of DMRG spin densities 

The dependence of the spin density distribution on the number of DMRG active-system 
states m is shown in Figure [5] where the differences in spin density distribution are plotted 
for DMRG(13,y)[m] calculations with respect to the converged DMRG(13,29)[2048] reference 
calculation. For increasing m-values, the differences in the spin density distribution decrease 
(see each row in Figure [5] from the left to the right). Similarly, we observe that the spin 
density gradually converges with respect to the dimension of the active orbital space (see 
last column from the top to the bottom of Figure [5j). In particular, changes in the spin 
density are negligible when m is increased from 1024 to 2048, and hence, reliable spin density 
distributions can be obtained even if the total energy is not yet converged with respect to 
m (the difference is 1.195 mH, see above). 
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Table III: Ground state energy for [Fe(NO)]^+ surrounded by four point charges at a distance 
of rfpc = 1.131 A from the iron center in Hartree atomic units for our DMRG(x,^)[m] calcu- 
lations employing different numbers of DMRG active-system states m. The CAS(11,14)SCF 
energy is —1393.013 396 Hartree. 

£'/Hartree 

Method DMRG(13,20) DMRG(13,24) DMRG(13,29)~ 

m = 128 -1393.014 662 -1392.991 085 -1393.014 010 

m = 256 -1393.018 626 -1393.019 309 -1393.024 883 

m = 512 -1393.020 065 -1393.021 876 -1393.030 374 

m = 1024 -1393.020 511 -1393.022 946 -1393.033 001 

m = 2048 -1393.020 646 -1393.023 294 -1393.034 196 



Furthermore, the Aabs- and Arg-values quantify the convergence series of the determined 
DMRG spin density distributions. In Table [IVl both error quantities are listed for each 
DMRG(x,^)[m] spin density with respect to the DMRG(13,29)[2048] reference spin density. 
In general, the absolute error Aabs and the root-square error A^s decrease for increasing 
m keeping the dimension of the active orbital space fixed. Note that larger active orbital 
spaces require a larger m-value to obtain the same accuracy as achieved in smaller active 
space calculations. This is not immediately evident from the error data presented in Ta- 
ble [IV] since diflFerent dimensions of the active orbital space are compared which result in 
nonzero error values, while error values determined for different parameter sets, but the 
same dimension of the active orbital space could vanish. The large error values for the 
DMRG(13,24)[128] calculation indicate that important states were not picked up by the 
DMRG algorithm resulting in the large differences in the spin density distribution displayed 
in Figure [5l Furthermore, since both error values determined for the DMRG(13,29)[1024] 
calculation are below the threshold values, no considerable improvement in the accuracy of 
the spin density distribution can be expected when m is further increased to more than 2048 
active-system states. 

In order to demonstrate the convergence of the DMRG(13,29) wave function with respect 
to m (and thus the convergence of the obtained DMRG(13,29)[2048] reference spin density 
distribution), the C ASCI- type wave function expansions are reconstructed and compared 
for all m- values. In particular, the infiuence of the missing rfa;2_^2 -double-shell orbital can 
be assessed by examining the CI coefficients corresponding to Slater determinants with an 
occupied rfa.2_^2 -double-shell orbital. Following the conclusions of a benchmark study for 
intermediate CAS sizes (see Supporting Information), only the most important configura- 
tions (I (7/ 1 > 0.00001) are necessary to obtain an accurate wave function expansion. As 
convergence threshold for the sampling procedure, a value of 0.001 is sufficient. With this 
threshold, similar CASCI-type wave function expansions are obtained for a quantum fidelity 
measure close to 0.998. The set of quantum fidelity measures F^.,mi+i for our five DMRG 
calculations with m, G {128, 256, 512, 1024, 2048} is {0.991800,0.995510,0.996983,0.997639}. 
As the number of DMRG active-system states is enlarged, the CI coefficients of the recon- 
structed wave function expansion converge gradually which is indicated by the increasing 
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Table IV: The absolute error Aabs and the root-square error Ars of the DMRG(13,2/)[m] 
spin densities with respect to the converged DMRG(13,29)[2048] reference spin density for 
a different number of normalized active-system states m for [FeNO]^^ surrounded by four 
point charges at a distance of rfpc = 1.131 A from the iron center. The Aabs and Ars values 
of the CAS(x,2/)SCF calculations of Ref.^° with respect to the DMRG(13,29)[2048] reference 
spin density are also listed. 



Method 


^abs 


Ars 


DMRG(13,20)[128] 


0.030642 


0.008660 


DMRG(13,20)[256] 


0.020088 


0.004930 


DMRG(13,20)[512] 


0.016415 


0.003564 


DMRG(13,20)[1024] 


0.015028 


0.003162 


DMRG(13,20)[2048] 


0.014528 


0.003028 


DMRG(13,24)[128] 


0.590022 


0.235922 


DMRG(13,24)[256] 


0.020993 


0.003245 


DMRG(13,24)[512] 


0.014045 


0.003633 


DMRG(13,24)[1024] 


0.011622 


0.002668 


DMRG(13,24)[2048] 


0.010731 


0.002361 


DMRG(13,29)[128] 


0.032171 


0.010677 


DMRG(13,29)[256] 


0.026005 


0.006790 


DMRG(13,29)[512] 


0.010826 


0.003406 


DMRG(13,29)[1024] 


0.003381 


0.000975 



CAS(11,11)SCF 
CAS(11,12)SCF 
CAS(11,13)SCF 
CAS(11,14)SCF 
CAS(13,13)SCF 
CAS(13,14)SCF 
CAS(13,15)SCF 
CAS(13,16)SCF 



0.086658 
0.080249 
0.046303 
0.042544 
0.052239 
0.073400 
0.053157 
0.104928 



0.024495 
0.020591 
0.011402 
0.010954 
0.012124 
0.019850 
0.011180 
0.031922 



quantum fidelity measure. Note that Fmi,m2 is close to the ideal value of 0.998 already for 
a small number of DMRG active-system states m, and hence, only minor variations in the 
large CI coefficients occur when m is increased which explains the slight differences in the 
spin density distributions displayed in Figure El^c) . 

To demonstrate that this is indeed the case, the CI coefficients of the most important Slater 
determinants {\Ci\ > 0.0001) corresponding to the m = 128 and m = 1024 calculations 
are shown in Figure IHl Slater determinants with large CI weights (|C/| > 0.05) are similar 
for both DMRG parameter sets, only minor deviations can be observed. Note that all these 
Slater determinants have been incorporated in the DMRG wave function already for m = 128. 
Considerable differences in CI weights are present for Slater determinants corresponding to 
small- valued CI coefficients (|C/| < 0.015), while some Slater determinants with |C/| < 0.01 
have not been incorporated in the DMRG wave function for m = 128. These off- size or 
missing configurations lead to the different spin density distributions for small-m values. 
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From the reconstructed CASCI-type wave function, the influence of the (ia;2_2^2-double-sheU 
orbital as weU as of the empty hgand orbitals on the spin density distribution can be analyzed. 
In the upper part of Table El conflgurations containing an occupied rfa;2_^2 -double-shell 
orbital and corresponding to the largest CI coeflicients are presented. In the lower part of 
Table El some selected conflgurations with large CI coeflicients carrying excitations to empty 
ligand orbitals that can not be included in standard CASSCF calculations are presented 
and compared for the DMRG(13,29)[128] and DMRG(13,29)[1024] calculations. In general, 
Slater determinants with an occupied rfa;2_^2-double-shell orbital feature small CI weights (< 
0.003) and are hence of minor importance, while Slater determinants bearing occupied ligand 
orbitals feature large CI coeflicients. Conflgurations containing occupied ligand orbitals that 
are only included in the DMRG(13,29)[m] calculations (marked in bold face in Table Ej) 
possess considerably large CI weights. All other Slater determinants with excitations to 
different empty ligand orbitals have smaller CI coeflicients. Hence, those ligand orbitals 
pose a signiflcant contribution in obtaining accurate spin density distributions for the small 
model complex and cannot be neglected from the active orbital space. 

Table V: Some important occupation number vectors (ONV) with the corresponding CI 
weights from DMRG(13,29)[m] calculations for [FeNO]^+ surrounded by four point charges 
at a distance of rfpc = 1.131 A from the iron center. Upper part: ONVs containing an 
occupied rfa;2_^2 -double-shell orbital (marked in bold face). Bottom part: additional selected 
important conflgurations with occupied natural orbitals that cannot be included in the active 
orbital space in CASSCF calculations (marked in bold face), for the same DMRG(13,29)[m] 
calculations. 2: doubly occupied natural orbital; a: natural orbital occupied by an a-electron; 
b: natural orbital occupied by a /3-electron; 0: empty natural orbital. 

CI weight 



Slater determinant m = 128 m = 1024 

b2b222a0a0000000 0000000 a 00000 0.003 252 0.003 991 

bb2222aa00000000 0000000 a 00000 -0.003 226 -0.003 611 

222220ab00000000 0000000 a 00000 -0.002 762 -0.003 328 

ba2222ab00000000 0000000 a 00000 0.002 573 0.003 022 

b2a222a0b0000000 0000000 a 00000 -0.002 487 -0.003 017 

202222ab00000000 0000000 a 00000 0.002 405 0.002 716 

b222a2a0b0000000 0000000 0000a 0.010 360 0.011 558 

22b2a2a0a0000000 0000000 bOOOO 0.009 849 0.011 366 

22b2a2a0b0000000 0000000 aOOOO -0.009 532 -0.011 457 

b2222aab00000000 0000000 0000a -0.009 490 -0.010 991 

a2222baa00000000 0000000 0000b -0.009 014 -0.010 017 

b2b222a0a0000000 0000000 OaOOO 0.008 820 0.010 327 

b2222aab00000000 OOaOOOO 00000 -0.004 277 -0.005 436 

22b2a2a0b0000000 aOOOOOO 00000 -0.004 224 -0.006 852 
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4.2 Assessment of CASSCF spin densities 

The converged DMRG(13,29)[2048] reference spin density can be used to assess the accu- 
racy of CASSCF spin density distributions and benchmark the quahty of the (restricted) 
active orbital spaces in standard CASSCF calculations (see Figure [7](a)). Note that the 
same isosurface value has been taken to display the DMRG(13,2/)[m]^DMRG(13,29)[2048] 
and CAS(x,2/)SCF-DMRG(13,29)[2048] spin density difference plots. The CASSCF spin 
density distributions determined for medium-sized active orbital spaces oscillate around the 
converged DMRG spin density. Depending on which double-d-shell orbital is included in the 
active orbital space, the /3-electron density around the NO ligand is either overestimated 
or underestimated. This results either in pure spin-polarized cases with /3-electron den- 
sity found only around the nitrosyl ligand for CAS(11,11), CAS(11,14), CAS(13,13), and 
CAS(13,14), or some additional a-electron density present around the NO ligand associ- 
ated with a simultaneous decrease in the /3-electron density for CAS(11,12), CAS(11,13), 
CAS(13,15), and CAS(13,16). 

Similarly, the large Aabs- and Ars-values stress the differences in the spin density distributions 
which are considerably larger than those from the DMRG(13,2/)[m]^DMRG(13,29)[2048] dif- 
ference analysis (Table HVl) . Furthermore, Table HVl indicates that the CAS(11,11)SCF- and 
CAS(11,12)SCF calculations and the CAS(11,13)SCF- and CAS(11,14)SCF calculations, re- 
spectively, are of similar accuracy, as they have similar error values, but the spin density 
difference plots emphasize the qualitatively different spin density distributions. Increasing 
the dimension of the active orbital space results in even larger deviations from the DMRG 
reference spin density because the active space is not stable and important orbitals are ro- 
tated out of the CAS. Note that aU DMRG calculations— except DMRG(13,24) [128]— yield 
smaller error- values and smaller differences in the spin density difference plots. 

Although the CASSCF spin densities are quantitatively converged with respect to the active 
orbital space, significant qualitative — but also non-negligible quantitative — differences to 
the DMRG(13,29)[2048] reference spin density can be observed. The extension of the active 
orbital space by including an additional shell of rf-orbitals only is not sufficient to obtain 
a qualitatively accurate spin density distribution for the small iron nitrosyl molecule. Our 
analysis indicates that empty ligand orbitals are essential for calculating reliable reference 
spin densities. This may have severe implications for the standard CASSCF approach that 
require further analysis in future work. 



4.3 Comparison to DFT spin densities 

A comparison of DFT and CASSCF spin density distributions for medium-sized active or- 
bital spaces for the [FeNO]^+ molecule has been already discussed in our previous work 
(see Ref."^ for more details). For an unambiguous benchmark of approximate exchange- 
correlation density functionals, the DFT spin densities of Ref.^^ can be compared to the 
DMRG reference distribution. The qualitative analysis of the DFT-DMRG(13,29)[2048] 
spin density difference distributions is shown in Figure [71(b) . When comparing to the results 
obtained in Ref . ,^21 similar conclusions concerning the performance of approximate exchange- 
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correlation density functionals can be drawn. The best agreement is found for BP86, BLYP, 
and TPSS, while the remaining approximate exchange-correlation density functionals yield 
larger deviations and result in too large spin polarization. We should note that BP86, BLYP, 
and TPSS correctly predict the distribution of the a-electron density around the nitrosyl 
ligand, although it is overemphasized. In general, nonhybrid functionals yield spin densities 
which are in closest agreement with the DMRG reference distributions. This observation is 
supported by both error measures which are smallest for BP86, BLYP, and TPSS (see Table 
I in the Supporting Information). 



5 Conclusions and Outlook 

In this work, we have demonstrated how reliable ab initio spin density distributions can be 
calculated for very large active spaces. Our procedure is based on the DMRG algorithm 
and on two different approaches to obtain spin density matrix elements: (i) on-the-fly di- 
rectly from the second-quantized DMRG elementary operators or (ii) from an approximate 
C ASCI- type wave function expansion which is determined by our SRC AS algorithm. ^^ The 
reconstructed CASCI-type wave function can also be used as a means to compare a series of 
DMRG calculations employing a different number of DMRG active-system states m. 

The small noninnocent molecule [FeNO]^+ surrounded by four point charges represents a 
suitable system to validate our approach. The spin density distributions are highly sensitive 
to the nature of the converged state. We deliberately converged DMRG wave functions that 
correspond to local minima in the electronic energy in order to compare with qualitatively 
wrong wave functions. The possibility of convergence into local minima is shown by exam- 
ining the (largest) CI coefficients of the SRCAS-reconstructed CASCI-type wave function. 
Strong deviations with respect to the absolute value of the CI coefficients indicate that the 
number of DMRG active-system states m is chosen too small, and hence important states 
have not been incorporated by the DMRG algorithm. Spin densities corresponding to such 
local minima deviate considerably from the ground state spin density. 

The convergence analysis of the spin density distribution for the [FeNO]^^ molecule consid- 
ered active orbital spaces comprising up to 29 active orbit als. Difference plots of the spin 
density distribution for different active orbital spaces as well as the absolute error and the 
root-square error in the spin density difference distribution indicate a quantitatively con- 
verged spin density with respect to the dimension of the active orbital space and the number 
of active-system states m (which was as large as tti = 2048). The DMRG reference spin den- 
sity has been used to validate CASSCF spin densities resulting in significant quantitative 
and even qualitative differences. Considering an additional shell of d-orbitals is not sufficient 
to obtain reliable spin densities for the small model system and the active orbital space must 
be extended by additional unoccupied ligand orbitals. Similar difficulties are likely to be 
present for larger iron nitrosyl complexes where the point charges are replaced by different 
ligands and hence additional ligand and iron orbitals must be included in the active orbital 
space. The DMRG study of larger {FeNO}^ complexes is now pursued in our laboratory. 

A convergence analysis of the spin density in terms of spin density difference plots with 
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respect to the number of DMRG active-system states indicates that rehable reference spin 
densities can be obtained even if total energies are not converged with respect to m. A 
similar conclusion was found in our previous work regarding the energy splittings of states 
of different spin multiplicity.'^'^'^ Comparison of CI weights corresponding to the most 
important configurations of the reconstructed CASCI-type wave functions for diflFerent m- 
values furthermore ensures that reliable spin densities are obtained. The similarities in 
DMRG wave functions can be quantified by the quantum fidelity measure which can be used 
as an additional convergence criterion for spin density distributions in a sequence of DMRG 
calculations. 

Spin densities calculated from approximate CASCI-type wave functions are in good agree- 
ment with the DMRG reference spin density. Qualitatively reliable spin densities can be 
obtained even for large thresholds of COM (0.001) when the most important configurations 
have been picked up in the wave function expansion. For this threshold, the CASCI-type 
wave function contains Slater determinants with CI weights larger than 0.00001 which are 
important for the spin density. 

The comparison of DFT spin densities with the DMRG reference distributions allows us 
to benchmark approximate exchange-correlation density functional. Although nonhybrid 
functional yield spin density distributions closest to the DMRG reference, significant qual- 
itative and quantitative differences to the DMRG reference distributions could be observed 
for all investigated density functional. Similar conclusions were drawn in our previous study, 
where DFT spin densities were assessed against CASSCF spin densities ,23 entailing that none 
of the investigated exchange-correlation density functional yields sufficiently accurate spin 
density distributions for the [FeNO]^+ molecule. 
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Figure 4: Mutual information and single orbital entropies 5(1) for the DMRG(13,y)[64] 
calculations determined for different numbers of active orbitals in the [Fe(NO)]^^ molecule 
surrounded by four point charges at a distance of rfpc = 1.131 A from the iron center. 
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Figure 5: DMRG(13,^)[m] and CAS(a;,y)SCF spin density difference plots with respect 
to the DMRG(13,29)[2048] spin density distribution (d) for [FeN0]2+ surrounded by four 
point charges at a distance of dpc = 1.131 A from the iron center. AU spin densities 
are displayed for an isosurface value of 0.001. (a) DMRG(13,20)[m]-DMRG(13,29)[2048] 
spin density difference plots, (b) DMRG(13,24)[m]-DMRG(13,29)[2048] spin density differ- 
ence plots, (c) DMRG(13,29)[m]-DMRG(13,29)[2048] spin density difference plots, (d) The 
DMRG(13,29)[2048] reference spin density distribution. 
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Figure 6: Distribution of the absolute value of the CI coefficients for the DMRG(13,29)[m] 
calculations with m = 128 and 1024, respectively, for [FeNO]^+ surrounded by four point 
charges at a distance of rfpc = 1.131 A from the iron center. The CI coefficients reconstructed 
for both DMRG calculations are always printed for the same Slater determinants. The 
determinants are ordered according to the CI weight of the DMRG(13,29)[2048] reference 
calculation. 
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Figure 7: (a) CAS(x,7/)SCF and (b) DFT spin density difference plots with respect to the 
DMRG(13,29)[2048] spin density distribution for [FeNO]^+ surrounded by four point charges 
at a distance of rfpc = 1.131 A from the iron center. All spin densities are displayed for an 
isosurface value of 0.001. 
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Figure 8: Distribution of the absolute value of the CI coefficients for the DMRG(7,7)[m] 
calculations employing different renormalized active-system states m and the CAS(7,7)SCF 
reference calculation for the [Fe(NO)]^"'" molecule surrounded by four point charges located 
at a distance of 1.131 A from the iron center. The CI coefficient obtained in the different 
calculations are always printed for the same Slater determinants. The determinants are 
ordered according to the CI weights of the CAS(7,7)SCF calculation. 
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Figure 9: CAS(7,7)CI[m]-DMRG(7,7)[m] spin density difference plots determined for differ- 
ent sampled subspaces of the complete A^-particle Hilbert space for the [Fe(NO)]^^ molecule 
surrounded by four point charges located at a distance of 1.131 A from the iron center. The 
spin density differences are plotted for the reconstructed CAS(7,7)CI-type wave function 
with respect to the DMRG(7,7)[m] reference spin density for the same number of renormal- 
ized active-system states m. The threshold (thr) indicates the upper bound of the COM 
value in the sampling-reconstruction procedure for each CASCI-type wave function. An 
isosurface value of 0.0003 was chosen for all spin density difference plots. 
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Figure 10: CI histogram of the absolute values of the CI coefficients for the Slater deter- 
minants for reconstructed CASCI-type wave function expansions from the DMRG(7,7)[7tz] 
calculations with different renormalized active-system states m for the [Fe(NO)]^^ molecule 
surrounded by four point charges at a distance of rfpc = 1.131 A from the iron atom. The 
CAS(7,7)SCF reference calculation is also shown for comparison, thr corresponds to the 
threshold value of COM in the sampling-reconstruction procedure and denotes the accuracy 
of the reconstructed CASCI-type wave function. All Slater determinants with CI coefficients 
in an interval as indicated on the abscissa are grouped together. 
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Figure 11: Excitation histogram for CAS(7,7)CI-type wave function expansion reconstructed 
for DMRG(7,7)[m] calculations with m renormalized active-system states and for different 
threshold (thr) values of COM for the [Fe(NCr)]^+ molecule surrounded by four point charges 
located at a distance of 1.131 A from the iron center. The CAS(7,7)SCF reference excitation 
pattern is also shown and marked as ref. 
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Figure 12: CAS(7,7)CI[m]-DMRG(7,7)[m] spin density difference plots determined for differ- 
ent sampled subspaces of the complete A^-particle Hilbert space for the [Fe(NO)]^+ molecule 
surrounded by four point charges located at a distance of 0.598 A from the iron center. The 
spin density differences are plotted for the reconstructed CAS(7,7)CI-type wave function 
with respect to the DMRG(7,7)[m] reference spin density for the same number of renormal- 
ized active-system states m. The threshold (thr) indicates the upper bound of the GOM 
value in the sampling reconstruction procedure for each GASGI-type wave function. An 
isosurface value of 0.0003 was chosen for all spin density difference plots. 



34 



0.1(5 


X 1 1 1 1 1 1 1 1 1 

- ^© — O Q^^^ 


1 1 


z 




G-Om=32 


0.01 








X m=48 


- 


0.001 




^ 0.0001 




% 


J 




_ 


y 




_ 


o 


I 






- 


• ^ 










.9 le-05 


— 






— 


^ 










«^ 










<D 


- 






- 


o le-06 


— 






— 


HH 


- 






- 


^ le-07 


_ 




_ 


<1> 


- 




- 


^4— > 










^ 


- 






- 


^ le-08 


— 






— 


C/3 










X) 


I 






I 


^ le-09 


— 






— 


le-10 


- 




- 


le-11 


- 




- 


le-12 


1 1 1 1 1 1 1 1 1 


1 1 (o 


~ 



1 10 

Slater determinants 

Figure 13: Slater determinants which are not incorporated by the DMRG algorithm for dif- 
ferent renormalized active-system states m in the DMRG(7,7)[7n] calculations for [Fe(NO)]^^ 
surrounded by four point charges in a distance of rfpc = 0.598 A from the iron atom. The 
missing Slater determinants are ordered according to their CI weight. Increasing m, embeds 
some of the previously missing configurations. The incorporated configurations for m = 48 
are displayed as open circles in the m = 32 curve, while crossed circles correspond to miss- 
ing Slater determinants for both m = 32 and m = 4:8. Note that important configurations 
which contain large CI weights could not be incorporated by the DMRG algorithm in small- ttz 
calculations. 
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CASSCF and DMRG spin density convergence for in- 
termediate CAS 

In Ref. [36], different medium-sized active orbital spaces are discussed ranging from 11 
electrons correlated in 9 orbitals up to 13 electrons correlated in 16 active orbitals. In our 
DMRG study, we first focus on five different active orbital spaces, namely on a CAS (11, 9), 
CAS(11,11), CAS(11,12), CAS(11,13), and a CAS(11,14). In the CAS(11,9), the active 
orbital space is composed of five Fe 3(i-orbitals {dxy^ dxz^ dy^^ rf^2, and dx2_y2)^ both NO 
TT- and 7r*-orbitals. The latter four CASSCF calculations contain an additional shell of 
Fe rf-orbitals: the Fe dxz- and d^y- (CAS(11,11)), dyz- (CAS(11,12)) and rf^2 -double-shell 
orbitals (CAS (11, 14) which contains also an additional empty ligand orbital). As orbital 
basis, natural orbitals from the corresponding CASSCF calculations are employed [68-70]. 
All CASSCF calculations as well as the calculation of the one-electron and two-electron 
integrals were performed with the Molpro program package [58] using Dunning's cc-pVTZ 
basis set for all atoms [59,60]. As the convergence of DMRG calculations can crucially depend 
on the ordering of orbitals in the algorithm [Moritz, Hess, Reiher, J. Chem. Phys. 2005, 
122^ 024107], the DMRG orbital orderings were optimized by minimizing the entanglement 
among molecular orbitals following the ideas of Refs. [49,57,72]. Additionally, we applied 
the CI-DEAS procedure in the start up procedure in order to keep the system size small 
and accelerate convergence [57]. All DMRG calculations are performed with the Budapest 
DMRG program [71]. The number of DMRG active-system states m ranges from 220 for the 
smallest active orbital space CAS (11, 9) up to 1680 for the CAS (11, 14) calculation in order 
to reproduce the corresponding CASSCF reference energy, while for the intermediate active 
orbital spaces m = 790 (CAS(11,11)), m = 960 (CAS(11,12)) and m = 1280 (CAS(11,13)) 
active-system states, respectively, are required to reproduce the CAS(x,^)SCF energies and 
spin density distributions. 

For all these DMRG calculations, the spin density distribution of the corresponding CASSCF 
reference calculation could be exactly reproduced and is therefore not displayed here. As 
an additional convergence criterion for the DMRG calculations, the reconstructed CASCI- 
type wave function of each DMRG calculation can be analyzed and compared to the ex- 
act CASSCF reference wave function. To perform such an analysis, we calculated the 
weighted quantum fidelity for a reconstructed CASCI-type wave function and the corre- 
sponding CASSCF reference wave function considering only the most important configura- 
tions (CI > 10~^). The quantum fidelity was further normalized with respect to the norm of 
the truncated CASSCF reference wave function since all configurations with CI coefficients 
smaller than 10~^ were neglected. The set of quantum fidelity measures Fdmrg,casscf for 
{CAS(11,9),CAS(11,11),CAS(11,12),CAS(11,13),CAS(11,14)} is {0.999999,0.999999,0.999999,0.999998,0.9 
Hence, the overlap of both wave functions is excellent, and similar spin density distributions 
are obtained. 
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Table VI: The absolute error Aabs and the root-square error A^s of DFT spin densities with 
respect to the converged DMRG(13,29)[2048] reference spin density for [FeNO]^^ surrounded 
by four point charges at a distance of rfpc = 1.131 A from the iron center. 



Method 


Aabs 


Ars 


OLYP 


0.360309 


0.0826634 


OPBE 


0.395280 


0.0891197 


BP86 


0.158527 


0.0346259 


BLYP 


0.138796 


0.0304860 


TPSS 


0.178942 


0.0425946 


TPSSh 


0.718517 


0.1634330 


M06-L 


0.446750 


0.1048820 


B3LYP 


0.951955 


0.2159410 
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